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pected single-telescope trigger rates would result in substantial data rates. To ameliorate these 
data rates, a novel, hardware-level Distributed Intelligent Array Trigger (DIAT) is envisioned. A 
copy of the DIAT operates autonomously at each telescope and uses reduced metadata from a 
limited subset of nearby telescopes to veto events prior to camera readout. We present the results 
of Monte-Carlo simulations that evaluate the efficacy of a "Parallax width" discriminator that can 
be used by the DIAT to efficiently distinguish between genuine gamma-ray initiated events and 
unwanted background events that are initiated by hadronic cosmic rays. 
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1. Introduction 

CTA |[T]] is a next-generation array of ground-based 7 -ray telescopes that will build upon the 
experience of current Imaging Atmospheric Cherenkov Telescope Arrays (lACTs) including VER¬ 
ITAS, H.E.S.S. and MAGIC. CTA has been designed to achieve unprecedented sensitivity, angular 
resolution and energy resolution over a broad range of 7 -ray energies. CTA will be a heteroge¬ 
nous array, comprising multiple telescope designs that will contribute variously to the overall per¬ 
formance characteristics of the array. A highly innovative telescope design considered for CTA 
uses a dual-mirror Schwarzschild Couder (SC) configuration, which incorporates a finely pixe- 
lafed camera comprising 11,328 independenf silicon phofomulfiplier-based readouf channels. The 
SC felescope (SCT) design provides a wide field of view (~ 8 °), while high-resolufion capfure of 
Cherenkov images provides excellenf angular and energy resolution. 

To ameliorate fhe effecl of Nighf Sky Background (NSB) photons randomly friggering each 
of fhe 11,328 individual readouf channels, and fhereby enable a reduction of fhe frigger fhreshold, 
lower resolution super-pixels (2,832) are formed from four adjacenf imaging pixels. The combined 
signal amplifudes of fhe super-pixels are used to form a Boolean-valued trigger image wifh summed 
oufpuf fhreshold level corresponding to 3.1 phofoelecfron-equivalenf (hereafter p.e.) counfs segre- 
gafing frue (friggered) and false values. An individual telescope is deemed to have generated a valid 
frigger condifion if ifs Boolean frigger image includes 3 or more adjacenf friggered super-pixels. 
lACTs operafe in a sfrongly background-dominated regime. Array friggering schemes use infor- 
mafion from mulfiple telescopes to vefo background evenfs before camera readouf wifh fhe goal of 
sfabilizing fhe array’s energy fhreshold and dead lime under variable ambienl illumination. Eor fhe 
densely pixelaled SCT camera, conlrol over fhe array frigger rate is also desirable fo guarantee Ihaf 
dala-lransfer rales remain Iraclable af fhe exlremes of normal observing conditions. 

Al low incidenl phofon intensifies, fhe overwhelming majorily of single- and super-pixel frig- 
gers are associaled wifh fhe ambienl NSB. The rate of random super-pixel friggers is often sufficienf 
to generale a large rale of spuriously valid individual felescope friggers. This lype of background 
can be efficienlly suppressed using a simple multiplicity array trigger lhal uses lemporal informa¬ 
tion fo refain only Ihose evenfs for which mulfiple telescopes frigger wilhin a nominal coincidence 
window. Eor larger photon inlensilies, anolher background componenl comprising lemporally cor- 
relafed Cherenkov lighl from cosmic-ray (CR) inifialed air showers becomes dominant. 

Electromagnetic air-showers initiated by 7 -rays are characterized by a single shower axis and 
triggered telescopes capture coherent elliptical images. In contrast, the hadronic component of 
cosmic-ray showers produces multiple sub-showers and typically yields a more fragmentary dis¬ 
tribution of Cherenkov light. The parallax width discriminator (P) leverages the difference 
between CR and 7 -ray shower images to rapidly distinguish between these event categories at the 
hardware level. 

2. Computation of the Parallax Width Discriminator 

It is instructive to separate the distinct algorithmic stages that comprise the computation of P 
into two separate categories. Eigure [I] illustrates the manner in which telescope-specific camera 
dafa are processed fo derive a compacf geomefrical represenfafion of fhe capfured shower image. 
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Figure 1: The left-hand panel illustrates an unprocessed, finely pixelated image corresponding to a 7 -ray- 
initiated air-shower. The central panel illustrates the intermediate, coarsely pixelated image that is derived 
by combining the signals from sets of 4 neighbouring imaging pixels. The right-hand panel illustrates the 
Boolean-valued trigger image that is required by the algorithm that computes P. The orange arrows corre¬ 
spond to the vector rp, which connects the implicitly unweighted trigger-image centroid with the fiducial 
camera-plane coordinate r* = (0,0) at the camera centre. For comparison, the white arrow {left-hand panel) 
connects the signal-amplitude-weighted mean position of all imaging pixels with r*. 
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Figure 2: Array mirror plane projections illustrating P computation for an on-axis 7 -ray event {top-left), an 
off-axis 7 -ray event with Ey<2 TeV {top-right), an off-axis 7 -ray event for Fy > 2 TeV {bottom-left), and 
an off-axis cosmic-ray proton event {bottom-right). 
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1. For each telescope that triggers, the centroid-vector rcj of the coarsely-sampled trigger image is 
defined as the mean camera-coordinates of the set of triggered super-pixels rcj = ((■rc,trig), (jc,trig))- 

2. For each centroid vector, a second vector rf.,- = re,; — r* is defined where r* is an arbitrarily 
selected fiducial camera-plane coordinate r* = (v:^,y^). To simplify computation, r* is defined 
to be the camera centre (0,0). 

The events shown in each panel of Figure ^ provide representative examples of y-ray and proton- 

initiated events with various simulated characteristics. They illustrate how the subsequent algorith¬ 
mic operations are used to compute P for each event class. 

1. The direction of each computed rp vector is projected from the independent coordinate systems 

of each telescope camera plane into a unified coordinate system spanning the mirror plane of 
the telescope array, defining a set of projected vectors The mirror plane is defined to 

intersect the mean geographical coordinate of each telescope comprising the complete array 
and to lie perpendicular to their common pointing axes. 

2. A set of nj mirror-plane coordinates {r/y} is computed that corresponds to the forward inter¬ 
sections between the projected vectors. 

3. Finally, the parallax width is defined as the dispersion of the intersection coordinates. 


P = 



( 2 . 1 ) 


The top-left, top-right and bottom-left panels of Figure ^ all correspond to y-ray-initiated events. 
For such events, the computed trigger-image centroids rc,/ correspond closely with the camera- 
plane projection of a 3-dimensional point on the air-shower axis at the height of maximal Cherenkov 
emission. The fiducial camera coordinate r* = (0,0) is also the camera-plane projection of a sec¬ 
ond three dimensional point. For y-ray-initiated events, there is a single, well defined shower axis 
and each rf.,- connects projected 3-dimensional points that are the effectively identical for all tele¬ 
scopes. Accordingly, the mapping {r^,} —)• yields a set of vectors which intersect at a tightly 

clustered region of the array mirror plane and the computed value of P is small. 

Figure ^ {top-left) corresponds to a y-ray-initiated events aligned with the telescopes’ optical 
axes (on-axis). When projected into the camera plane’s coordinate system, the arrival directions of 
on-axis events correspond precisely with the fiducial camera-plane coordinate r* = (0,0), which 
is implicitly the camera-plane projection of a 3-dimensional point on the air-shower axis. Con¬ 
sequently, the derived {rj j} coincide closely with the projection of the axis into the array mirror 
plane. The top-right and bottom-left panels represent ojf-axis y-ray events with energies below and 
above 2 TeV, respectively. For lower energy events, the 3-dimensional-point pairs connected by 
each of the remain identical for all telescopes and the tight clustering of {rj j } is preserved. At 
higher energies, only one of the projected points maps to the air shower axis, so the coincidence 
between the {f/y} and the mirror-plane projection of the air-shower axis is lost. 

Finally, the bottom-right panel illustrates a proton-initiated event. There is now no guarantee 
that Yc4 for each telescope corresponds with the projection of a single 3-dimensional point, since 
different telescopes may image multiple, different sub-showers of the hadronic cascade. Accord¬ 
ingly, the tight clustering of the {rj y} is lost and the computed value of P is large. 
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E = 1.9662 TeV ^ 16 Trigger Pixels E = 1.9662 TeV ^ 16 Trigger Pixels 




Figure 3: Computation of the super-pixel multiplicity nxp that is required to satisfy the pass-through cri¬ 
terion. In the left-hand panel, the energy-binned effective collection area for profon-initiated events that 
trigger at least one telescope is folded with the expected cosmic ray spectral shape (oc and used to 

determine the typical proton energy Epyo% at which the expected rate of single telescope triggers falls below 
10% of its peak value^. In the right-hand panel, the energy-binned distribution of trigger image super-pixel 
multiplicities is used to derive the value of nxp that corresponds to events for which Ep ~ Ep iq%. 

3. Pass-Through Trigger for High-Energy Events 

The parallax width algorithm assumes that the centroids of trigger images reliably encode 
the geometry of the air-showers. If the ground-plane intensity profile of Cherenkov light that is 
emitted by y-ray-initiated air-showers is substantially asymmetric, the validity of this assumption 
may degrade. Ideally, the trigger-image centroid rc,,- should correspond closely with the. full image 
centroid r^ defined as fhe signal-amplifude-weighfed mean position of all imaging pixels fhaf 
frigger in response fo incidenf Cherenkov phofons. Wifhouf access fo fhe information provided by 
fhe individual pixel amplitudes, fhe Boolean-valued trigger images appear more symmefric and fhe 
rc,i thaf are used fo compufe P may nof accurafely represenf fhe air-shower geomefry. 

The lefi and right-hand panels of Figure [I] illusfrafe schemafically how any camera-coordinafe 
offsefs Ec.i = rc,i ~ i befween fhe fwo cenfroid definifions produce corresponding direcfional 
perfurbafions of each rf./- Mapping fhese misaligned vectors yields a sef fhaf fypically in¬ 

creases fhe dispersion befween fhe {r/y} infersecfion coordinafes, and inflafes fhe compufed value 
of P. The pofenfial magnifude of ecp increases for high-energy y-ray-inifiafed air showers, which 
fypically produce exfensive images fhaf comprise a large number of friggered super-pixels. 

To prevenf spurious rejecfion of genuine high-energy, y-ray-inifiafed evenfs, a pre-calibrafed 
mulfiplicify fhreshold utp is used fo unconditionally accepf (or pass through) evenfs for which any 
felescope in fhe array capfures a frigger image comprising utp or more super-pixels. As illusfraled 
by Figure fhe expecfed single-telescope trigger rate Rp for simulated, proton-initiated air-showers 
is used to calibrate an appropriate value for njp- To retain effective suppression of the most frequent 
cosmic-ray triggers, a threshold corresponding to the typical super-pixel multiplicity for proton 
events that trigger at 10% of the peak single-telescope rate is adopted. 

^Note that y-ray-like images of a hadronic sub-shower typically sample a fraction of the energy of the incident 
proton. Accordingly, genuine y-rays that produce images comprising «xp trigger pixels have energies that are typically 
~ 30% of protons that do so. 
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Figure 4: Top row: The computed distributions (left) and CDFs (right) of P values corresponding to point- 
origin (on-axis, red) and diffuse (blue) 7 -ray events, and diffuse cosmic-ray background events (green). 
Bottom row: Collection areas for a point-origin 7 -ray source assuming nxp —^ °° (left), and htp = 16 (right). 

4. Results and Discussion 

To investigate the efficacy of P to distinguish 7 -ray and cosmic-ray initiated air-showers, simu¬ 
lations of the SCT array configuration indicated by the magenta markers in Figure § were produced 
using the sim_telarray software package^ |^. The raw simulated camera images were pro¬ 
vided as input to a computer simulation of the DIAT array triggering scheme outlined above. Figure 
^ (top-left) shows the three distributions of P values that are computed for each investigated dataset, 
while the right-hand panel displays the corresponding cumulative distributions for P. For the point¬ 
like 7 -ray source, all the simulated photons are incident on-axis, and 90% of computed P values 
are < 11 m, in accordance with expectation. For extended 7 -ray sources the majority of incident 
photons are incident off-axis. Nonetheless, the expected prevalence of low P values is realized with 
90% of events that do not fulfil the pass-through criterion having P <21 m. In contrast ~ 75% of 
all proton-initiated events yield F > 27 m, confirming the expectation that P is a highly efficient 
discriminator between 7 -ray- and proton-initiated events. 

The bottom-left-hand panel of Figure ^ plots the energy-binned effective collection areas for 
the simulated array configuration, assuming a point-like 7 -ray source, without application of a 
pass-through for events yielding extensive high-multiplicity trigger images. Curves corresponding 
to four distinct array-triggering criteria are shown for comparison. A fiducial scenario, correspond- 

^ Simulations of point-like, extended (2° radius) 7 -ray emission, and a proton background were used. The 7 -ray 
events model astrophysical sources located centrally in the field of view at an altitude of 70° and an azimuth of 180°. 
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ing to a single telescope trigger requirement (i.e. no array trigger) is shown in black. The green 
curve illustrates the effective area for events that trigger at least two telescopes and are subsequently 
found to yield two high-quality Cherenkov images^. The violet curve illustrates the collection area 
for a more desirable subset of events that yield at least three high quality Cherenkov images. Com¬ 
putation of P requires at least two mirror-plane intersection coordinates {rj y}, which implicitly 
requires a telescope trigger multiplicity of 3 or more, but does not require that any of the triggered 
telescopes yield high-quality Cherenkov images. The red curve illustrates the collection area af¬ 
ter application of the parallax width criterion, rejecting any events for which P > 25 m. At low 
energies, the parallax-width and three-telescope trigger areas are comparable, but at higher ener¬ 
gies, the parallax-width trigger retains a fraction of minimally reconstructible events that yield less 
than three high-quality Cherenkov images. The bottom-right-hand panel of Figure Q illustrates 
how adoption of a pass-through threshold njp > 16 further enhances the retention of minimally 
reconstructible y-ray events by the DIAT. 

The top-left-hand panel of Figure ^ uses simulated effective collection areas to demonstrate 
the response of the DIAT to an extended astrophysical y-ray source, adopting njp > 16. At low 
energies the rejection of events yielding F > 25 m retains a large fraction of events that yield 
three high quality Cherenkov images, while at energies > 1 TeV, the pass-through trigger results in 
retention of all minimally reconstructible events. The top-right-hand panel of Figure ^ illustrates 
the power of P to effectively reject cosmic-ray proton-initiated events. The blue curve illustrates 
the effective collection area for a traditional two-telescope multiplicity array trigger, that accepts 
events for which at least two neighbouring telescopes trigger. Both array triggering strategies 
reduce the energy-integrated collection area by a factor ~ 7, but it is evident that the parallax- 
width trigger outperforms the two-telescope multiplicity requirement at lower energies where the 
incident cosmic ray rate is largest. Figure ^ {bottom row) illustrates how DIAT using the parallax 
width algorithm can be used to improve the low-energy sensitivity of the SCT subarray. Field- 
to-field variability of the NSB intensity between observations can induce unpredictable spikes in 
the array trigger rate if accidental coincidences between spurious single telescope triggers cannot 
be effectively suppressed. Typically, the required suppression is achieved by requiring a higher 
single-pixel trigger threshold, which increases the overall low-energy threshold of the array. The 
ability of the DIAT to perform real-time multi-telescope event vetoing allows lower pixel thresholds 
to be used while maintaining tractable array trigger rates. In Figure the differential trigger 
rates for a Crab-pulsar-like y-ray source are compared for events that satisfy P < 25 m assuming 
the nominal trigger pixel threshold of 3.1 p.e., with the set of events that fulfil a traditional two- 
telescope multiplicity array trigger. The lower pixel threshold made feasible by a hardware-level 
DIAT yields a substantial improvement in sensitivity below 200 GeV, achieving a factor of ~ 20 
enhancement at the lowest energies. 

In summary, simulations show that efficient rejection of cosmic-ray events is possible with 
a hardware array trigger using the moments of Cherenkov light images. The algorithm success¬ 
fully discriminates between background and genuine y-ray triggers, retaining a large majority of 

^High-quality gamma-ray images are those that are subsequently useable for parameterization and geo¬ 
metric reconstruction of the corresponding air shower and the properties of its progenitor y-ray. Hereafter, 
events yielding at least two high-quality images will be described as minimally reconstructible 
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Figure 5: Top row: Collection areas for diffuse 7 -ray emission (left), and diffuse protons (right) correspond¬ 
ing to a pass-through super-pixel multiplicity threshold njp = 16. Bottom row: The expected array trigger 
rates for on-axis 7 -ray-initiated events (left) are shown for a Crab-pulsar-like spectral shape. The red curve 
corresponds to application of the parallax width algorithm with P < 25 m, a pass-through threshold nxp >16, 
and a pixel trigger threshold corresponding to 3.1 p.e. per trigger pixel. For comparison, the blue curve plots 
the expected array rate for a traditional two-telescope multiplicity array trigger, with a trigger-pixel threshold 
of 4.5 p.e. Furthermore, the ratio of the red and blue curves (right) are shown, indicating that application of 
the parallax width discriminator enables effective recovery of low energy events. 


reconstructible events. By vetoing spurious single and multiple telescope triggers before data are 
read out, the algorithm reduces the array trigger rate, enabling finely sampled events with heavy 
data payloads to be generated by participating SCTs without overwhelming the array data trans¬ 
fer infrastructure. Real-time consideration of data from multiple telescopes also allows the rate 
of NSB-induced array triggers to be controlled without increasing the single-pixel trigger thresh¬ 
olds. The resultant enhancement of low-energy sensitivity may be particularly useful for studies of 
spectrally soft targets like the Crab pulsar. 
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